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CI Abstract The present work is a brief review of the progressive search of improper <5-functions 

which are of interest in quantum mechanics and in the problem of motion in General Rela- 
' tivity Theory. 

' A great deal of my work was just playing with equations and seeing what they give. 

1 P. A.M. Dirac 

+1» , Every one who works, no matter how briefly or superficially, in complex relativity 

will find himself acknowledging Jerzy's work. 
£S ' M.P. Ryan 

; 1 Introduction 

> • 

j^. | The advent of quantum mechanics opened a new domain of concepts, including generalized 

functions. Introduced as an improper function by P.A.M. Dirac in 1926 [1] (see also the 
Dirac's book [2]), the delta function had been used in physics for quite time before the 
formal work of L. Schwartz, published in 1950 [3]. The mathematical foundations of gen- 
eralized functions, however, appear to have first been formulated in 1936 by S.L. Sobolev 
in his studies on the Cauchy problem for hyperbolic equations [4] (see also the works of J. 
Hadamard [5] and M. Riesz [6]). At the present time, the distribution theory has advanced 
substantially and has found a number of applications in physics and mathematics. Indeed, 
the use of generalized functions leads to remarkable simplifications in the problems that 
one usually faces in contemporary physics. 

One of important applications of the Dirac's delta, out of the quantum area, occurred 
in general relativity. In 1927, A. Einstein and J. Grommer reported the first solution of 
the problem of motion; in the procedure, they used the delta function to represent matter 
in the field equations [7]. Thus, the simplification of the problem was done at the cost of 
introducing singular structures. Primary refinements were done in 1938 by Einstein, L. 
Infeld, and B. Hoffman [8]. Some years later, in their 1960 book Motion and Relativity, 
Infeld and J. Plebahski discussed a number of deeper improvements to the Einstein, Infeld, 
and Hoffman approach [9]. In particular, they modified the definition of the delta function 
in order to properly manage singularities in general relativity. 

In this paper we examine the origin of generalized -improper- functions in quantum 
and general relativity theories. Although history leads to a better understanding of the 



modern physics, it is not my interest to cover the entire development but rather to fix 
attention on some specific points. The paper is organized as follows: Section 2 overviews 
the origin of the delta function in quantum mechanics. Section 3 deals with the aspects of 
the motion problem in general relativity (GR) connected with the modifications developed 
by Infeld and Plebahski on the Dirac's delta function. Section 4 is devoted to some other 
"delta objects" appearing quite recently in the literature and to the concluding remarks. 

2 The Dirac's delta function 

One of the most interesting developments of quantum mechanics concerns the concept of 
commutativity and starts in 1925, with an idea of W. Heisenberg: one ought to ignore the 
problem of electron orbits inside the atom, and treat the frequencies and amplitudes asso- 
ciated with the line intensities as perfectly good substitutes. In any case, these magnitudes 
could be observed directly [10]. Indeed, it is necessary to bear in mind that in quantum 
theory it has not been possible to associate the electron with a point in space, considered as 
a function of time, by means of observable quantities. However, even in quantum theory it 
is possible to ascribe to an electron the emission radiation [11] (see Van der Waerden [12] 
p 263). He was certain that no concept enter a theory which has not been experimentally 
verified at least to the same degree of accuracy as the experiments to be explained by the 
theory [13]. Thereby, Heisenberg concluded that the physical variables should be repre- 
sented by specific arrays of numbers (matrices). A conclusion which, in turn, led him to 
an apparently unexpected result: Whereas in classical theory x{t)y{t) is always equal to 
y{t)x{t), this is not necessarily the case in quantum theory [11] (see Van der Waerden [12] 
p 266). It was almost inconceivable that the product of physical quantities could depend 
on the multiplication order. 

Thus, in exchange for the classical notion of position and momentum in atoms, Heisen- 
berg introduced the concept of observables (measurable experimental magnitudes) and re- 
marked on their non-commutation properties. After the approval of W. Pauli, Heisenberg 
published his results in the paper Quantum theoretical reinterpretation on kinematic and 
mechanical relations (Uber quantentheoretische Undeutung Kinematischen und mechan- 
ishen Beziehunge) [11] and gave a copy to M. Born. When Born read the paper he noticed 
that the Heisenberg's symbolic multiplication was the matrix algebra. Later on, Jor- 
dan, Heisenberg, and Born published the three men's paper On Quantum Mechanics (Zur 
Quantenmechanik) [14] in which they reported a matrix formulation of the new theory 
(See also Van der Waerden [12], Jammer [15], Duck-Sudarshan [16], and Mehra [17,18]). 

On 28 July 1925, during a stay in Cambridge with R.H. Fowler, Heisenberg delivered 
the talk "Term Zoology and Zeeman Botany" before the Kapitza Club. The subject 
dealt with the anomalous Zeeman effect and the enormous difficulties to build atomic 
spectroscopy up by means of ad hoc rules, a remarkable topic for somebody who had solved 
the quantum puzzle recently (see Mehra [18] Ch 19.10). A month later, Heisenberg sent 
the proof-sheets of his paper to Fowler who, in turn, gave it to his research student, Dirac. 
After reading the paper, Dirac pondered it for two weeks and noticed that Heisenberg's 
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idea had provided the key to the 'whole mistery' [18]. In his own words: non- commutation 
was really the dominant characteristic of Heisenberg's new theory [19]. Dirac concluded 
that quantum mechanics could be inferred from the Hamilton's form of classical dynamics 
by considering new 'canonical variables' obeying a non-commutative 'quantum algebra'. 
The results were published by Dirac between 1925 and 1927. One of his papers, The 
physical interpretation of the quantum dynamics [1], was decisive in the formalization of 
the new theory. There, an arbitrary function of the position and momentum is shown 
to be smeared over the entire momentum space if the position is infinitely sharp (the 
uncertainty principle!). The main point of the formulation was a transformation theory 
which required the introduction of the improper function 5(C)- 

Dirac reconsidered the Heisenberg's idea of observables: When we make an observation 
we measure some dynamical variable. . . the result of such a measurement must always be 
a real number. . . so we should expect a real dynamical variable. . . Not every real dynamical 
variable can be measured, however. A further restriction is needed (see Dirac [2] p 34- 
35). Then, he formalized the concept by defining an observable as a real dynamical 
variable whose eigenstates form a complete set, and stated that, at least theoretically, every 
observable can be measured. If the eigenvalues of the observable ( consist of all numbers 
in a certain range, then the arbitrary eigenkets \X) and \Y) of Q may be expressible as 
the integrals 



and \C") being eigenkets of ( belonging to the eigenvalues £' and Q" respectively, x and 
y labelling the two integrands, and the range of integration being the range of eigenvalues. 
We say that \X) and |y) are in the representation of the basic bras \0 (a similar definition 
is true for discrete eigenvalues). By considering the product (X\Y), we take the single 
integral 



The integrand in (|2j) vanishes over the whole range of integration except at the point 
C = Following Dirac's formulation (see Dirac [2] Ch 10), as in general (^|^) does not 
vanish, so in general {C' x \C"y) must be infinitely great in such a way as to make ([2]) non- 
vanishing and finite. To get a precise notation for dealing with these infinite objects, Dirac 
introduced the quantity 5(C), depending on the parameter £ and fulfilling the conditions 
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And, a more general expression: 
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The range of integration in Q and (JI|) does not need to be from — oo to oo, but may 
be over any domain 0, surrounding the critical point at which the 5 function does not 
vanish. Dirac acknowledged there is something unusual about the delta function and 
decided to call it an 'improper function'. The following expressions are essentially rules of 
manipulation for 5 functions 
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The first of equations ([5]) means that C^(C)i as a factor in an integrand, is equivalent 
to zero. The second one is easy to verify by using integration by parts n times, and 
means that 5(C) can be formally differentiated as many times as one wishes. There are 
diverse ways to face the 5-function. For example, it appears whenever one differentiates 
a discontinuous function like the Heavside onej 
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I i C>o, 

for which one gets dO(C) /dC, = 5(C)- Of special importance is the Fourier representation of 
<5(C)- It is obtained trough the eigenfunctions of the operator id/d(, that is (2ir)~ 1 ^' 2 e~ lk ^ , 
henceforth 



6(C'-C") = — [e ik «'-^dk. 
2vr J 
u 



(7) 



Thus, the 5-function is just a shorthand notation for limiting process which simplifies 
calculations. In general, one can take a class of functions 6(e, Q, such that 



limy /(cH(e,CK = /(o). 



(8) 



In practice, one uses the following mathematical scheme: all calculations have to be per- 
formed not on (5(C) but on S(e, £)• The limiting procedure e — > has to be made in the 
very last result. Two plausible models are the sequence of Gaussian distribution functions 



S(e,0 



1 



and the (simplest) set of square well potentials 

V(v,C) = \ "I 



-C 2 /2e 2 



ld<e 
ICI>e 
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(10) 



1 Indeed, 6(C) is also an improper function. Hence, the derivative dO/d(^ should be understood as in 
the context of equation 
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with v = — 1/(2 e); e > 0. Finally, the integrand of equation ([2]) will now be written 



were we have dropped the labels x and y. If we are interested in two different representa- 
tions for the same dynamical system P, the quantities (i]\Q are called the transformation 
functions from the representation {[??)} into {|C}}- The ket \P) will now have the two 
representatives {rj\P) and (Cl-P)j defining the corresponding transformation equations (see 
Dirac [2], p 75). According to Dirac, the transformation functions are example of prob- 
ability amplitudes. Thus, the statistical interpretation of Born is also applicable in the 
Dirac's transformation theory. 

Nowadays, the formulae ([3])-([8]) are easy to analyze but in Heisenberg-Dirac times 
neither matrix formalism nor improper functions were popular among theoreticians. As 
we have said, Heisenberg was advised by Born on the connection between matrix algebra 
and his non-commutative operations. In contradistinction, the Dirac's training as engineer 
seems to be fundamental; in his own words: all electrical engineers are familiar with the 
idea of a pulse, and the 5-function is just a way of expressing a pulse mathematically (see 
Jammer [15] p 316). However, although the 5- function is attributed to Dirac, it had been 
introduced in physics much earlier. Prior to its appearance in quantum mechanics, it 
was used by Hertz in statistical mechanics in connection with the concept of temperature 
(see Jammer [15] pp 317). Its occurrence in pure mathematics was noticed in 1815 by 
A.L. Cauchy whose derivation of the Fourier-integral theorem is based on the modern 
(5-function [20] (the derivation is reproduced in Van der Pol and Bremmer [21] Ch 8). 
Also in 1815, S.D. Poisson worked on the Fourier-integral theorem and followed a similar 
procedure as that of Cauchy [22]. In 1882, G.R. Kirchhoff used the Green's theorem in 
the study of Huygen's principle. Kirchhoff too was acquainted with the improper function 
delta, which he denoted by F: As to the function F we assume that it vanishes for all 
finite positive and negative values of its argument, but that it is positive for such values 
when infinitely small and in such a way that 



where the integration extends from a finite negative to a finite positive limit [23] . Kirchhoff 
remarked on the fact that 2<5(/i~ 1 /\/2, £) approximates F for very large \i (see eq. (J2])). In 
1891, influenced by the works of Cauchy and Poisson, H. Hermite proposed the integral [24] 



and analyzed its limit A — > for small values of 8. It is easy to see that by taking 8 = 0, 
provided a and (5 lie at different sides of t = 0, the above integral defines the class of 
functions S(X,t) converging to 5(t) in the sense of ([8]). Later on, the delta function is put 
forward in the work of O. Heavside [25] (see equation ([6]) and the Heavside books [26]). 

Apart of the previously quoted antecedents, the intuitive Dirac's procedure opened an 
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intriguing problem in pure mathematics. Some criticism was presented, e.g. von Neumann 
wrote: the Dirac's method adheres to the fiction that each self-adjoint operator can be put 
in diagonal form. In the case of those operators for which this is not actually the case, 
this requires the introduction of 'improper' functions with self-contradictory properties. 
The insertion of such mathematical 'fiction' is frequently necessary in Dirac's approach... 
It should be emphasized that the correct structure need not consist in a mathematical 
refinement and explanation of the Dirac method, but rather that it requires a procedure 
differing from the very beginning, namely, the relevance on the Hilbert theory of operators 
[27]. The first attempts to mathematically formalize the definition of the 5-function 
were done in 1926-1927 by Hilbert and published later by Hilbert, von Neumann and 
Nordhein under the title The foundations of quantum mechanics ( Uber die Grundlagen 
der Quantenmechanik [28].) The challenge was finally faced by Schwartz in the context of 
his theory of distributions [3] (see also Sobolev [4]). Thus, the ill-defined 5-function and 
its derivatives were replaced by well-defined linear functionals (distributions) which have 
always other distributions as derivatives on the test functions space. 

The Dirac's function plays an alternative role in quantum mechanics: it is an exactly 
solvable potential enjoying many useful applications [29]. As a physical model, it has 
been used to represent localized matter distributions or potentials whose energy scale is 
high and whose spatial extension is smaller than other relevant scales of the problem. 
Arrays of 5-function potentials have been used to illustrate Bloch's theorem in solid state 
physics (the Dirac's comb) and also in optics where wave propagation in a periodic medium 
resembles the dynamics of an electron in a crystal lattice. The bound state problem in 
one dimension (for potentials involving either attractive or repulsive delta terms) has an 
exact implicit solution whenever the eigenvalue problem without deltas can be solved 
exactly [30]. In two and higher dimensions it provides a pedagogical introduction to the 
techniques of regularization in quantum field theory [31] (in one dimension, the quantum 
system needs no regularization.) It has been also studied in the context of supersymmetric 
quantum mechanics [32-36] where the susy partner of the attractive (5-function is the 
purely repulsive (5-function [32,33,35]. Similar results are obtained for potentials made 
up of additive 5-function terms [34, 36] . 

3 The Infeld-Plebahski's delta function 

In all descriptions of nature we use two alternative concepts: field and matter. Matter is 
composed of particles and the field is created by moving particles. The picture is simple 
at the cost of having singular fields. Furthermore, what can one say about the motion of 
the sources? 

In Newtonian gravitational theory the concept of gravitational field is reduced to the 
action at a distance. However, according to relativity, no linear field theory can determine 
the motion of its sources because no action can be propagated with a speed greater than 
the speed of light. Hence, one must add the motion equations to the field equations. 
The statement is no longer true in nonlinear field theories as Einstein and Grommer have 



6 



shown in their paper of 1927, General Relativity Theory and Laws of Motion (Allgemeine 
Relativitatstheorie und Bewegungsgesetz) [7]. They obtained an unexpected conclusion: 
the equations of motion for a test particle are but a consequence of the field equations! 
The Einstein and Grommer paper opened as well more general problems in GRT. One of 
them was how to find whether the equations of motion of two particles can also be deduced 
from the field equations: a challenge which remained open till 1936, when L. Infeld arrived 
in Princeton to begin a collaboration with Einstein. 

Before his presence in Princeton, Infeld was working with M. Born in Cambridge. They 
faced the problem of modify Maxwell's electrodynamics so that the self energy of the point 
charge is finite. Their results are nowadays known as the Born-Infeld electrodynamics. As 
Einstein rejected the idea from the very beginning and Born and Infeld did not succeed 
in their attempts to reconcile it with quantum theory, Born 'warmly recommended' Infeld 
to Einstein (see Born [37] p 121) and, in fact, Infeld became Einstein's collaborator and 
assistant. Later on Einstein wrote to Born: We [Infeld and Einstein] have done a very 
fine thing together. Problem of astronomical movement with treatment of celestial bodies 
as singularities of the field (Born [37] Lett. 71, p 130). The 'fine thing' concerned a 
fundamental simplification of the foundations of GR. 

The Infeld-Einstein collaboration on the problem of motion was more persistent than 
any problem Infeld had tackled before: For three years I worked on this problem whose 
only practical application that I know of is the analysis of the motion of double stars by 
methods giving deeper insight than the old Newtonian mechanics. For three years I have 
been bothering with double stars [38]. In principle, the movement of sources is determined 
by the geodesic lines of the space-time world; the metrics of which satisfy the Einstein's 
field equations. The point of departure was the Einstein's assertion that the first part of 
this assumption is redundant; it follows from the field equations by going to the limit of 
infinitely thin, mass-covered world lines, on which the field becomes singular (see Born [37] 
p 131, and Infeld [39]). Infeld remarks the calculations were so troublesome that we decided 
to leave on reference at the Institute of Advanced Study in Princeton a whole manuscript 
of calculations for other to use [9]. 

Such a quantity of work was finally rewarded. In 1938 Einstein, Infeld and Hofffman 
(EIH) published the paper The Gravitational Equations and the Problem of Motion [8], 
in which the two-body problem was solved for the first time. However, as the relativ- 
ity non-linear field equations are too cumbersome to be solved exactly, approximation 
methods were required. The basic idea behind the EIH method is to take into account 
that for a function depending on coordinates and time, developed in a power series in the 
parameter c _1 , the time derivatives will be of a higher order than the space derivatives. 
In general, by using singularities to represent matter, the method consisted in forming 
certain two-dimensional surface integrals over surfaces enclosing these singularities. The 
field equations prescribed the laws by which the surfaces enclosing the singularities, and 
hence these singularities moved [39]. 

The Einstein-Infeld collaboration continued some more years bringing a progress in the 
problem of motion whose final solution will never influence our daily lives and will never 
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have any technical application. It is a purely abstract problem. An even more skeptical 
thought of Infeld was / do not believe that there are more than ten people in the world who 
have studied our papers on the problem of motion [38]. As it seems nowadays, Infeld had 
underestimated the importance of their own results. 

Not long after the EIH work was successfully completed; Infeld, then at the University 
of Toronto, published a paper with Wallace in which the EIH approach is applied to the 
problem of motion in electrodynamics [40]. Ten years later, in 1951, Infeld and Scheidegger 
worked on the problem of gravitational radiation reaction in the EIH formalism^). That 
same year Infeld left Canada [38, 42-44] and returned to Poland to join the Institute 
of Theoretical Physics in Warsaw University. Once in Warsaw, his work attracted the 
attention of an amount of brilliant graduate students and researchers all interested in 
gravitational wave theory. Bialynicki-Birula, Suffczyhski, Trautman, Werle, Plebahski 
and Krolikowski -the last two, students of Rubinowicz- are some of the names of the 
Infeld 's group. 

Almost fifteen years had elapsed since the EIH paper was published and the problem of 
motion attracted again the interest of Infeld. He worked on the subject with his group for 
about six more years to collect finally all their results in the book Motion and Relativity, 
written together with Jerzy Plebahski in 1960 [9]. The collaboration was profitable for 
young Plebahski; Infeld introduced him to non-linear electrodynamics (see for instance 
[45]), unitary operators [46] and to the problem of motion in GR among other topics. The 
elaboration of Motion and Relativity took almost four years of discussion on the contents 
and typeset. Both Infeld and Plebahski, finished the first chapter and appendiceal when, 
in 1957, Plebahski received a Rockefeller Fellowship to go to the United States. Before 
leaving, Plebahski prepared a sketch in Polish of the rest of the manuscript, except the 
last chapter which was added by Infeld lone [9, 42] (see also Garcfa-Compean et al [49] in 
this volume). 

The birth of the Infeld-Plebahski delta functions occurred at this phase of the research. 
In some sense, it followed the early Einstein's ideas: the energy— momentum tensor Tq>@ 
in the field equations introduced indeed an excess of information mixing the physics and 
geometry. The use of G a p = R a p — \g a p = —8-7rT al 3, with T a p proportional to the 
Dirac's function <W (the presence of matter), permits to skip the redundance, reducing 
the geometry to the singular solution of G a p = R a p- ^Saff R = [39]. Such interpretation, 
adopted by Infeld [50], permitted to simplify the entire deduction of the equations of 
motion. The concrete mathematical model was obtained in collaboration with Plebahski 
[9,51]. 

To fix the ideas, let £ s (t) be a world line and (p a scalar field that depends on coordinates 

2 Their paper Radiation and Gravitational Equations of Motion [41] gave rise to a considerable flow of 
discussion concerning the radiation reaction problem in GR, a subject in which Infeld was involved as soon 
as he arrived in Princeton. See the interesting Kennefick's paper [42] for details. 

3 Plebariski was almost 23 when Infeld arrived in Poland in 1950. Not long after they began in a close 
collaboration which strongly influenced the Plebahski's first research contributions. Indeed, it was in 1957 
that Infeld invited Plebahski to write their monograph. 

4 The main results of their appendix 1 "The S function", were published in 1956 in the paper On 
modified Dirac 8 -functions [47]. Later improvements were published in 1957 [48]. 
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x s , time x° = t, and also on the £, s (t) and their time derivatives £ s (i): <p = ip(x s ,t;^ s ,^ s ). 
As the procedure produces fields ip which are singular on the world lines £ s (i), Infeld and 
Plebahski looked for a transformation theory changing ip into a continuous function <p of 
the £ s , £ s and £ s , without recurse to the renormalization procedure. Hence, they were 
faced with the problem of interpreting the expression 

f 6 $dC, k>0 (12) 
n 

often considered as divergent. However, the diverse definitions of the ^-functions, as 
presented in the preceding section^, are useless to interpret (fT2|) . In all cases the integrands 
/(C) <KC) were considered for the continuous functions /(C), at least in the vicinities of 
C = 0. Thus, (I12j) is simply meaningless! Infeld and Plebahski solved the problem by 
narrowing the definition of Dirac's 5-function so that the integral 

/ HC)§ip(C)dC (13) 
h 

acquires a definite meaning even if ifj(Q has a singularity up to the kth order [47,48]. 
They introduced their delta function 5ip in an axiomatic form which extends the limits 
of Schwarts distribution theory: 

(IPi) Sip has all derivatives for C 7^ 0. 
(IP 2 ) Sip = if C ¥= 0; 5 IP (0) is undefined. 
(IP3) For a continuous function /(C) : 

/(0) = J f(()5ip(()d(. (14) 
h 

(IP4) For a certain k we have 

/^«="» < 15 » 

n 

where u>k is a previously assigned value. 

The axioms (IPi)-(IP4), all can hold for a convenient class of functions 6i P (e, C) in the 
realistic approach [9]. Hence, Infeld and Plebahski defined the following algebraic rules 

<5 /P (c) = (i + ^|Cl fc )<MC) 

5 IP (t) = a\t\ k ±(C5(C)) 

J Infeld and Plebanski categorize three different methods: (A) axiomatic, essentially depicted by prop- 
erties (J3j) — dSJ) ; (B) Fourier transformation. This lies on the definition (0 and the preceding properties; 
(C) The realistic method, which lies on the very core of the Dirac's intuition (see the paragraphs between 
equations j?} and ((9|) for which the axiomatic properties ([3])-([5]) are immediately justified [9,47,48]. 
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where 5jp corresponds to the choice cok = in (I15p . 5 is the conventional Dirac function 
and a is an infinite constant chosen such thati 



J 5 IP (()d( = l = J S IP (OdC 



The meaning of equations (|16p is that their two sides give equivalent results as factors in 
an integrand, exactly like the ordinary Dirac's ^-function. The (5/p-functions, however, 
allow to associate definite meanings with integrals of products of Sip with functions that 
become divergent for £ — > 0. This is why Infeld and Plebahski called them "good functions" 
[9]. They claimed the application of these functions as equivalent to the regularization 
procedure. 

As for the transformation theory, Infeld and Plebahski generalized these good functions 
to more dimensions (which is simple enough) and established [39] 

<P<*JP(3) (X S ~€ S )d(3)X = % 

as the definition of <p, where 5/p(3) is their three-dimensional good J-function. Thus the 
tilde means two things: singularities of tp are ignored and, for x s , the £ s are substituted. 

4 Is the 5— Zoology exhausted? 

It might seem that the <5-Zoology (by paraphrasing Heisenberg [17]) is exhausted. Yet, 
from time to time, distributions on differential domains are also considered. For in- 
stance, a rigorous definition of the delta function could be obtained, in the sense of 
Mikusihski [52,53], by defining the generalized functions as the closure of certain ordi- 
nary functional spaces with respect to a weak topology. Last years, the hyper functions 
of Sato [54] (considered more general than the improper functions) have more and more 
applications. In the present section, we shall analyze a set of new objects which have 
been profitable in the context of Darboux transformations in quantum mechanics (a topic 
which, under the name of factorization method, was also investigated by Dirac, Infeld, 
and Plebahski). Let us start by remarking that functions ([9]) and (|10p behave as e _1 for 
small values of e. Now, what about functions 8{e, Q showing an arbitrary e-dependence 
instead of Unlike the Dirac's case, we shall take a family of "well potentials" 

f -e" 2 ICI < e 

S(e,C := - (17) 

{ |C| >£• 

In order to get a wider meaning of this new objects, let us analyze 

e- 2 (e(-£,o) 

A'(-.u { -e- 2 Ce (Q,s) (18) 
|C|>e. 



6 This is certainly possible in the realistic method and a turns out to be as singular as £ fe . In their 
book, Infeld and Plebahski found a = e' h ^2 (ft+1)/2 r (^Mr 1 )] for a Gaussian S IP (e, (). 
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which has been reported by the Christiansen's group [55] and by Boykin [56] independently. 
It is a matter of integration (in the sense of distributions) to verify that lim £ _>o A'(e, Q = 
S'(C)- Thus A' approaches the derivative of the Dirac's delta function! The Christiansen's 
group worked on definition (jl8[) to analyze the scattering properties by regularizing finite- 
range potentials (point or contact interactions). Their approach leads to the conclusion 
that 5'(C) is a transparent potential as opposite to the Seba's theorem [57] which estab- 
lishes that <5'(C) should have zero transparency. They also studied second and third order 
derivatives of the delta function. On the other hand, Boykin obtained ()18f) by conveniently 
transforming a finite difference formula. He used a three dimension version of A' to derive 
the Gauss' law in a dielectric medium directly from the charge densities, without using 
potentials [56]. 

Now, as we can see, our well potentials (fT7|) resemble those in (fTHj) . In some sense, 
we could interpret either the limit 5(e, C)e^o as an "incomplete" derivative of the delta 
function or A'(e, £) as a combination of 5(e, £): 

A / (e,C) = 5(e,C)e(C)-<5(e,C)e(-C). 

Remark that, for dealing with test functions /, the sequence 5(e,Q guarantees small 
domains of integration f2, centered at the origin^. Henceforth, consider the following 
result 

J /(c) icfco *: = -/(£) (19) 

h 

where ^ £ (— e, e), the function /(£) is differentiable enough and the mean value theorem 
for integration has been applied. By calculating the limit e — > and interchanging the 
limiting process with integration in (|19p , one establishes the following rule of manipulating 
the (^-function: 

HO ICI = -s(C) (20) 

with 5 the ordinary Dirac's function. Now, let us draw our attention to the first of 
equations ([5]). It shows that whenever one divides both sides of an equation by a variable 
C which can take on the value zero, one should add on to one side an arbitrary multiple of 
5(Q) (see Dirac [2]), just as it occurs for the derivative of the log(C) function: dlog(() / d( = 
1/C — iir5((). Thus, equation (f20l) becomes: 

HC) = - 6 j^ + c5(C)sgn(C) (21) 

where c is an arbitrary constant and sgn(£) is the sign improper function. The right hand 
side of (|2ip is neither a definition of 5(C) nor an equality sensu stricto. As before, equations 
(I20p and (|21|) are merely operational equivalences requiring the integration. Nonetheless, 
5 has a stronger divergence than 5 at £ = 0. 

7 Indeed, f2 does not need to be centered at the origin but at an arbitrary accumulation point £o. I 
prefer to use £o = 0. 
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What can we accomplish with these new 5 'improper functions'? Let us take a 
continuous function (j>(C) which is not necessarily differentiable at £ = but such that 
(j)(C) ~ const\^\ for £ — > 0. Then, the following transformation holds: 

J mmd( = m (22) 

h 

where </>(£) is a new continuous and differentiable function in all the real line (remember 
we have taken the accumulation point Co equal to zero). As an immediate example one 
can substitute 4>(C) for f(() \(\ and 0(£) for — /(£) in equation (fT9|) ; after the usual limit 
procedure one gets ([22]) . 

Observe that, in the previous derivations, the class of functions 5(e, C) have been con- 
sidered on a free particle background. The relevance of the singular 5 'function' is analogous 
on a nontrivial background. In particular, let V(C) be a singular, one dimensional, po- 
tential growing as £~ 2 for £ — > 0. The new potential T4- eff (e,C) = V(C) + a^(e 5 C) can be 
proved to be regular at £ = for the appropriate value of the strength a and any value 
of e ^ 0. Recent results show that periodic singular potentials admit a regularization 
procedure in this sense [58]. If the initial potential is the Scarf's one: V S {Q =Vq/ sin 2 (£), 
then V* eg (s, C) is a family of regular potentials such that lim £ ^o Keg( e > C) = ^ S (C) + <KC)- 
Similar results are obtained for other singular potentials defined on the complete real line, 
including the cases of discrete spectrum (see e.g. Dutt et al [36], and Negro et al [59]). 
Furthermore, it has been shown that this procedure does not change the results of a super- 
symmetric transformation. Thus, the Darboux transformations and the <5-regularization 
procedure commute in quantum mechanics [58,59]. 

Finally, let us remark on the fact that generalized functions can be represented as 
sequences of ordinary functions which converge in a certain way. This property, as we have 
seen, stimulated the development of the theory of distributions and related approaches. 
The use of improper functions thus expands the range of problems that can be tackled in 
mathematical and theoretical physics, in particular in the theory of differential equations 
and quantum physics. 
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